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1. Introduction 



I The Cosmic Microwave Background Radiation (CMBR) provides a strong 

observational foundation for the standard cosmological scenario, the Big 
Bang theory. It is difficult to understand how to produce a 2.7°K black- 
lO ' body spectrum except in the context of the Big Bang scenario. The near 

. blackbody spectrum of the CMBR along with it's near isotropy provides 

I compelling evidence for a period of fairly quiescent Friedman-Robertson- 

' Walker expansion for many expansion time before recombination. The past 

Q^i decade has seen huge advances in the measurement of the CMBR, with 

' I COBE's definitive discovery of anisotropics and measurement of a near per- 

5^ ■ feet blackbody spectrum. The small deviations from isotropy have and will 

^ i continue to tell us a great deal about the inhomogeneities in our universe, 

' and small deviations from a blackbody spectrum can also tell us about the 

energetics in our universe. Such deviations have already been discovered in 
\ the direction of clusters of galaxies, although the mean CMBR spectrum 

^ • is, so far, indistinguishable from a blackbody spectrum. 

^ Here we give a introduction to the observed spectrum of the CMBR and 

discuss what can be learned about it. Particular attention will be given to 
how Compton scattering can distort the spectrum of the CMBR. This is 
left toward the end though. Unfortunately the author has no expertise in 
the area of how these measurements are made but Smoot has covered this 
area in his lectures. An incomplete bibliography of relevant papers is also 
provided. Some old but still highly useful reviews of the physics behind the 
spectra are by Danese and DeZotti[20], and Sunyaev and Zerdovich[74]. 
Theoretically not much has changed in this field in over 25 years. Much of 
the interesting work was done by Zel'dovich and Sunyaev in 1969[83]. 
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2. Executive Summary 

The universe today is fairly diffuse and cold, however the universe is ob- 
served to be expanding, and in the past we may deduce that the uni- 
verse was more dense and because of p dV work the matter in the universe 
would also have been hotter. Extrapolating the expansion back to very early 
epochs the universe would have been very hot and very dense and the uni- 
verse must have been expanding very rapidly in order to have grown as large 
as it is observed to be. Hence the Hot Big Bang. When the matter in the 
universe is hot and dense the thermal cqTiilibration time becomes very short. 
Thus we expect everything to rapidly approach thermal equilibrium and we 
therefore expect the photons in the universe to have a thermal (blackbody) 
spectrum at early times. It is easily shown that expansion of the universe 
(or traversal through any gravitational field) leaves a blackbody spectrum 
a blackbody spectrum although the temperature may change. This tem- 
perature change is known as the redshift and can sometimes be thought 
of as either a Doppler shift or a gravitational redshift. Formally speaking 
the two may be thought of as the same phenomena and there is often no 
physical sense in trying to separate them. 

Thus as a first approximation we expect the photons in the universe 
to have a blackbody spectrum. The fact that the cosmic microwave back- 
ground radiation (CMBR) has nearly a blackbody spectrum is strong evi- 
dence for the Hot Big Bang hypothesis. There is simply not enough matter 
around today to thermalize so many photons (there are > 10^ photons for 
every atom) and in any case most of the matter in our universe is much 
hotter that 3 K. The reason we might expect a deviation from blackbody is 
because some of the matter in the universe has gone out of thermal equilib- 
rium with the photons and may either heat or cool the photons. This can 
be done by non-equilibrium scattering or absorption of existing photons or 
by non-equilibrium emission of new photons. Clearly most of the matter 
in the universe is not today in thermal equilibrium with the CMBR and 
the spectrum offers us a probe of this. However there are so many more 
CMBR photons in the universe than there are protons or electrons that it is 
difficult for the matter to significantly distort the spectrum of the CMBR. 
Thus the fact that the observed CMBR spectrum is so close to a blackbody 
should come as no surprise. 

3. Measures of Temperature 

The brightness or specific intensity of light, I^, is defined as the incident 
energy per unit area, per unit solid angle, per unit frequency, per unit time. 
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It may be written 

= —^n^ (1) 

where v is the frequency and n-y(z^) is the quantum-mechanical occupation 
number, i.e. the number of photons (in each polarization state) per unit 
phase space volume measured in units of . Here h is Planck's constant, 
and we assume the light is not (linearly) polarized so that there are an 
equal number of photons in each polarization state. A blackbody or Planck 
spectrum has 

< = ^1^^ (2) 



exp 



where T is the the temperature. The high-frequency (hi' 3> kT) limit of 
the Planck spectrum is known as Wien 's law: 

C = ?^expf-^) (3) 



c2 kT 

while the low frequency (hv <^ kT) limit of the Planck spectrum is known 
as the Rayleigh-Jeans law: 

= . (4) 

Note that the intensity is proportional to the temperature in this case. 
One may invert the Planck spectrum and characterize the intensity by the 
thermodynamic temperature or brightness temperature: 

(5) 



kin 



2/11/3 

1 + 



Occasionally radio astronomers may define the brightness temperature by 
it's Rayleigh Jeans limit: 

In the radio region this is an excellent approximation to the thermodynamic 
temperature and is simply related to the intensity, and is therefore closer 
to what is actually measured. 

Here we are interested in small deviations from a blackbody spectrum, 
i.e. we have some temperature Tj which is a good fit to Tb at many frequen- 
cies, and want to express the actual spectrum in terms of small deviations 
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from a blackbody with this temperature, in particular in terms of the devi- 
ations in intensity from the blackbody spectrum, A/,/. For small deviations 
the deviation in brightness temperature is 

ATb = Tb - = ^ , / -r^ ^ = TTTT ■ 7 

Experimentally it is often easier to measure differences rather than abso- 
lute numbers: differences in intensity in different directions on the sky, or 
between the sky and internal calibrators. Note that in the Wien region dif- 
ferences in brightness temperature are greater than in the Rayleigh-Jeans 
region for the same difference in intensity. In what follows we will tend to 
plot spectral distortions in terms of differences in brightness temperature 
versus the dimensionless frequency, x. These "derived" quantities are prob- 
ably more appealing to a theorist than an observer since they are further 
removed from what is actually measured. 



4. Measured Mean Spectrum of the CMBR 

Over the years there have been many measurements of the CMBR. There 
have been many claims that the spectrum deviated significantly from a 
blackbody, especially in the Wien region, however recent measurements 
with FIRAS (Far-InfraRcd Absolute Spectrophotometer) on the COBE 
satellite has shown definitively that the spectrum is very close to a black- 
body [48]. Contemporary and quite competitive with with the first FIRAS 
measurements was a rocket experiment [34] which also found a blackbody. 
The most recent FIRAS results have appeared in ref [30] which are plotted 
in figs 1 &; 2. In the 2nd figure the we have converted to a more theoreti- 
cal representation by plotting versus x = and converting to fractional 
changes in temperature. The reason that this transformation is useful is 
that we can predict the shape of the deviations from blackbody in terms of 
X, which we cannot in terms of f since we have no a priori knowledge of T^. 
To transform to an x variable one must decide on a fiducial temperature. 
We have used the best-fit temperature, = 2.728K, taken from the most 
recent results of FIRAS[30]: = 2.728±0.002K. Note that the uncertainty 
in is which is larger than the error bars on most of the individual points 
in the plots. The reason that the uncertainty in Tj, — T-y can be smaller 
than the uncertainty in Ty is that the experiment measures the difference 
in brightness between a blackbody and the sky. The zb0.002K represents 
the uncertainty in the temperature of the reference blackbody, while the 
accuracy to which this reference is thought to be a blackbody is much bet- 
ter than this. Note that when fitting for a distortion from a blackbody one 
must fit for both the amplitude of the distortion and for simultaneously. 
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Figure 1. Plotted are the residuals in Rayleigh- Jeans temperature from the best fit 
blackbody as a function of frequency as stated by Fixsen et al. (1996). The error bars 
are 1-cr. The solid line is the subtracted Galaxy model at the Galactic poles. We see 
that these measurements are running up against a fundamental limitation of Galactic 
contamination. 



While FIRAS certainly revolutionized the field, and does make obsolete 
most other short wavelength measurements of the CMBR spectrum, it only 
looked at the frequency range 68 — 640GHz. The bolometric techniques used 
by FIRAS only work at high frequencies and therefore the spectrum at low 
frequency was not touched by FIRAS. There has been ongoing measure- 
ments of the absolute CMBR flux in the Rayleigh-Jeans region for 30 years 
and we list some results from the last 15 years in Table 1. One can see that 
measurement did not stop after COBE. As we shall see some of the spec- 
tral distortions we are looking for are most visible in the Rayleigh-Jeans 
regime. We have selected some of the most sensitive of measurements to 
plot in fig 3. The uncertainties vary widely with frequency and are orders 
of magnitude larger than those of FIRAS. Several authors have noted that 
these low frequency measurements tend to indicate a temperature lower 
than that obtained at higher frequencies [62, 7], suggesting that there may 
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Figure 2. Same as fig 1 except here we plot the fi-actional deviation in brightness 
temperature vs. the dimensionless frequency x = ^ . To do this we have used the best-fit 
photon temperature T-y = 2.728 K. Plotting things in this way accentuates the deviations 
at high frequencies. 



be a deviation from a blackbody spectrum at low frequencies. 

Measurement of the absolute CMBR spectrum, at the present level 
of sensitivity, face significant problems of Galactic contamination at both 
long and short wavelengths. Synchrotron radiation contaminates the long- 
wavelength spectrum while the short wavelength region is contaminated by 
dust emission. Since we cannot expect to observe the CMBR from outside 
of the Galaxy this is a fundamental limitation. Many of the results plotted 
here include significant corrections for this contamination. While there is 
a limit to how well one can subtract off the Galaxy, we can look forward 
to improvements in Galaxy modeling using results from anisotropy exper- 
iments which will have increasingly better sensitivity, sky coverage, and 
angular resolution. While anisotropy experiments cannot generally make 
absolute measurements of intensity, they can help to map out the Galaxy. 
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TABLE 1. Listed are measurements, made over the past 15 years, of the absolute 
CMBR brightness temperature at a variety of wavelengths. The results often 
include significant corrections for Galactic emission. Millimeter wavelengths are 
omitted as they have been superseded by results of FIRAS (> 68 GHz - see 
Fixsen et al. 1996). The ADS code given refers to the paper where these results 
are presented or reviewed and may be used to find the papers and abstracts online 
at the NASA Astrophysics Data System and mirror sites: adsabs.harvard.edu, 
cdsads.u-strasbg.fr, d01.ratk.nao.ac.jp . These codes are of the form year journal 
volume page. 



Frequency Wavelength Tcmbr, 1st Author ADS Bibliographic 

(GHz) (cm) (Kelvin) Code 



1.47 


20.4 


2.26 


+0.19 
-0.19 


Bensadoun 


1993ApJ. 


..409.. ..IB 


90. 


0.22 


2.60 


+0.09 
-0.09 


Bersanelli 


1989ApJ. 


..339..632B 


2.0 


15. 


2.55 


+0.14 
-0.14 




1994ApJ. 


..424..517B 


3.7 


8.1 


2.59 


+0.13 
-0.13 


De Amici 


1988ApJ. 


..329..556D 


3.8 


7.9 


2.64 


+0.07 
-0.07 




1990ApJ. 


..359..219D 


3.8 


7.9 


2.64 


+0.06 
-0.06 


?) 


1991ApJ. 


..381..341D 


25. 


1.2 


2.783 


,+0.02,5 
'-0.025 


Johnson 


1987ApJ. 


..313L...1J 


7.5 


4.0 


2.60 


+0.07 
-0.07 


Kogut 


1990ApJ. 


..355..102K 


1.410 


21.26 


2.11 


+0.38 
-0.38 


Levin 


1988ApJ. 


..334... 14L 


7.5 


4. 


2.64 


+0.06 
-0.06 




1992ApJ. 


..396....3L 


4.75 


6.3 


2.70 


+0.07 
-0.07 


Mandolesi 


1986ApJ. 


..310..561M 


2.5 


12. 


2.79 


+0.15 
-0.15 


Sironi 


1986ApJ. 


..311..418S 


0.600 


50. 


3.0 


+ 1.2 
-1.2 


5) 


1990ApJ. 


..357..301S 


2.5 


12. 


2.50 


+0.34 
-0.34 


51 


1991ApJ. 


..378..550S 


0.82 


36.6 


2.7 


+ 1.6 
-1.6 


?7 






2.5 


12.0 


2.78 


+0.3 
-0.3 


Smoot 


1987ApJ. 


..317L..45S 


33.0 


0.909 


2.81 


+0.2 
-0.2 


?i 




■>■> 


1.41 


21.2 


2.22 


+0.55 
-0.55 


V 




7) 


3.66 


8.2 


2.59 


+0.14 
-0.14 


75 




55 


10. 


3.0 


2.61 


+0.06 
-0.06 


55 




75 


90. 


0.33 


2.60 


+0.10 
-0.10 


55 




?5 


1.4 


21. 


2.65 


+0.33 


Staggs 


IQQ.SPhDT fiS 


-0.30 






10.7 


2.80 


2.73C 


1+0.014 
-0.014 


55 


1996ApJ. 


..473L...1S 


90. 


0.33 


2.57 


+0.12 
-0.12 


Witebsky 


1986ApJ. 


..310..145W 



5. Spectral Distortions of the CMBR 



While one should not be surprised that the CMBR has close to a blackbody 
spectrum, there are various mechanisms which should cause deviations from 
a thermal spectrum. Now we discuss a few of them. 
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Figure 3. Plotted are a selection of the low frequency measurements of the CMBR 
brightness temperature listed in table 1. Prom left to right the points are from Sironi 

et al. (1990), Bcnsadoun et al. (1993), BcrsancUi et al. (1994), Sironi& BonclU (1986), 
DeAmici et al. (1991), Mandolcsi et al. (1986), Kogut et al. (1990), Levin et al. (1992), 
Smoot et al. (1987), Staggs (1996), Johnson & Wilkinson (1987) and wore chosen because 
of the small orrorbars. The black band at the right indicates the FIRAS data (Fixson 
et al. 1996), while the horizontal straight line represents a temperature 2.728 K given by 
the FIRAS best fit blackbody spectrum. The long-dashed represents a chemical potential 
distortions with amplitude u = ±9 x 10~^ while the solid line gives free-free distortions 
with amplitude Yg = ±10~ . These are both idealized curves and one may expect (model 
dependent) corrections long- ward of 10 GHz (see Burigana, DeZotti, and Danese 1995). 



5.1. ANISOTROPIES 

The most common way in which the CMBR spectral distortion occurs is 
when the photons have a blackbody spectrum in each direction but the tem- 
peratures characterizing these spectra arc different in different directions. 
This direction dependent temperature difference is called anisotropy. The 
anisotropy can either be caused by Doppler/gravitational effect or because 
the gas emitting the photons really did have different temperatures. 

The first anisotropy discovered was the dipole anisotropy, i.e. the tem- 
perature varies like the cosine of the angle from some point on the celestial 
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Figure 4- Superimposed on the FIRAS data of fig 2 is the largest chemical potential 
distortion allowed by the data (Fixsen et al. 1996): fi = ±9 x 10~^. The falling positive 
curve is the far more plausible positive chemical potential distortion and the negative 
rising curve is a negative chemical potential distortion. 



sphere. It is usually attributed to the Sun moving at 371km/s. Note that 
to a first approximation, when averaging over the sky, the dipole does not 
contribute to the mean spectrum of the CMBR. 



One way to check that measured anisotropics are what they are sup- 
posed to be is to measure the spectrum. For a small anisotropy the change 
in flux from the mean spectrum should be proportional to the derivative of 
the flux of a blackbody with respect to temperature. FIRAS has done just 
that for the dipole[29, 30] and found just what was expected. Most modern 
anisotropy experiments use many frequency channels in order to check the 
spectrum of the anisotropy, or more specifically to be able to subtract off 
contamination of the measurements by other effects than the anisotropy. 
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Figure 5. Superimposed on the FIRAS data of fig 2 is the largest j/-distortion allowed 
by the data (Fixsen et al. 1996): y = ±1.5 x 10~^. The rising curve is the more plausible 
positive j/-distortion and the falling curve is a negative y-distortion. 



5.2. CHEMICAL POTENTIAL DISTORTIONS 

There are three processes which are important from thermahzing the CMBR 
spectrum in the early universe: Compton scattering, double Compton scat- 
tering, and free-free scattering (also known as bremsstrahlung) . Compton 
scattering is a much more rapid process but since it conserves the number 
of photons so it can only thermalize the energy distribution of the pho- 
tons and not the number of photons. All of these processes become more 
efficient as one goes to earlier and earlier epochs and eventually photon 
non-conserving processes start to become important. 

There is a epoch between z = 10^ and z = 2x 10^ during which Compton 
scattering is efficient in thermahzing the energy distribution while other 
processes are not capable of thermahzing the photon number. During this 
epoch, if the energy-to-photon ratio is perturbed from that required for a 
blackbody spectrum, the spectrum will instead approach a Bose-Einstein 



11 



distribution 

= (J \ 

where and ji ( the dimensionless chemical potential) are determined by 
the total energy available and the total number of photons available. If one 
starts out with a thermal distribution of photons at temperature Tj and 
injects a fractional increase in the energy density, without significantly 
increasing the number of photons one obtains 











0.456 - - 


) 







AU 

— = 0.714// <C1 (9) 



Since one must fit the observations to both T-y and /x one really can only 
measures the - term. Double Compton scattering and free-free scattering 
become increasingly more efficient at lower frequencies and there are usu- 
ally corrections to this formula at small frequencies x <^ 1[9, 10]. These 
corrections are not liable to be important for FIRAS measurements. 

This distortion to the spectrum is greatest at small x, however the 
FIRAS measurements at high frequencies are so accurate that they yield 
much better constraints on /x than does the low frequency experiments. 
Comparing with the FIRAS data one finds < 9 X 10-5 at the 2a level[30]. 
We compare the maximal allowed distortion to the low &i high frequency 
data in figs 3&;4, respectively. 

Thus we find the extremely stringent constraint at a very early epoch 

^ < 6 X 10-5 10^ < z < 2 X 10^ (10) 

Of course this is not to say that one expects large energy injection at these 
epochs. 

Note that for z > 10^ the CMBR spectrum is not telling us much about 
the energetics of the universe. However one can use Big Bang Nucleosyn- 
thesis to probe the total energy of the universe up to z ~ lO^'^ . 



5.3. Y DISTORTIONS 

If energy is injected into the universe after z ~ 10'^ Compton scattering 
is unable to thermalize the distribution. The fact that we observe very 
little deviation from a blackbody spectrum tells us that not much energy 
could have been injected compared with the thermal energy of the CMBR. 
If a small amount of energy is injected then one may solve for the linear 
perturbation from a blackbody spectrum under the action of Compton 
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scattering as was done by Zel'dovich and Sunyaev[83]. One finds that the 
perturbation in the photon occupation number is 



iVe is the number density of free-electrons, and ctt is the Thomson cross- 
section. If one could manage to cool gas below the radiation temperature 
one could produce a distortion with negative y, but typically this distortion 
is produced by ionized gas which is much hotter than the photons. In the 
early universe when the density of electrons is large even a small heating of 
the gas over the photon temperature may lead to a significant distortion. 
This distortion is generally referred to as a y-distortion when applied to 
the mean CMBR spectrum, but is usually called the Sunyaev-Zel'dovich 
distortion when referring to an anisotropy in the spectrum because there 
is more or less hot gas in one direction than another. Large amounts of 
hot ionized gas exist in clusters of galaxies and the "S-Z effect" has been 
observed in the directions of several clusters. There is no evidence for a y of 
the mean CMBR spectrum, although with sensitive enough measurements 
we should see the hot gas we know is out there. Fixsen et al.[30] have 
placed a limit of \y\ < 1.5 x 10~^ from the FIRAS data. We compare the 
maximal distortions to the FIRAS data in fig 5. Note that a positive y- 
distortion produces a negative change in at low frequencies and positive 
change in T^, at high frequencies, just what one expect if one was heating a 
fixed number of photons. This negative AT^ is sometimes called the "S-Z 
decrement" , for the S-Z effect was first looked for at radio wavelengths. 

5.4. DISTORTION FROM FREE-FREE 

Another important process for the CMBR spectrum is free-free scattering, 
which is the scattering of a free electron off of a charged nucleon either 
emitting or absorbing a photon; in most cases of interest emitting rather 
than absorbing. This is the same process which produces the X-rays ob- 
served from hot cluster gas operating at microwave and radio frequencies. 
The effect of free-free scattering on the CMBR spectrum is mostly likely to 
be seen at the very longest wavelengths measured. In this Rayleigh-Jeans 
limit the distortion it produces may approximated by[3] 




(11) 



where the "y-parameter" is 




(12) 





(13) 
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Figure 6. Plotted is the constraint on the temperature of a fully ionized universe as a 
function of redshift. The horizontally hatched region is excluded since \y\ < 1.5 x 10~^ 

while the diagonally hatched region would be excluded if Yg < lO"**. The sohd hne 
indicates where the gas temperature equals the photon temperature and the dashed hne 
gives the temperature as a function of redshift for a model where the gas is ionized by 
very massive stars (VMOs - see Stebbins & Silk 1986). The cosmological parameters used 
are Ho = 65km/sec/Mpc, f2o = 0.4 and fib = 0.10. 



where 



K = 



SmcikT^)^V6TTmekTe 



(14) 



and the Gaunt factor, ^ ~ 2 in most cases of interest. Note the Te ^ factor 
in K which means that the effect is suppressed for higher temperature gas. 
The 1/x^ dependence on this distortion means that it is the low frequency 
measurements which wih constrain it's amphtudc. In fig 3 we plot free-free 
distortions with Yg = ±10^^. We see that this size distortion is close to 
what is being constrained by these measurements, although no proper sta- 
tistical analysis has been done. This size free-free distortion would produce 
no significant effect in the FIRAS data, although if one went far enough 
into the Wien tail one would find large distortions from free-free emission. 
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Figure 7. The same as the previous figure except with f2o = 1 and f2b = 0.06. The 
constraints are less severe for for larger f2o and smaller f2b. 



In figs 6&7 wc have used constraints on y and Yq to put constraints 
on the temperature and epoch of a reionized universe, assuming presently 
favored cosmological parameters. We see that it really the y-distortion 
which is most important, the free-distortion only being detectable on the 
off chance that the gas was ionized and cold. Even though the limits on y 
are quite small we see that there is not too much of a constraint of ioniza- 
tion after z ~ 100. The constraints could be made stronger if one assumes 
a larger baryon density or a smaller total density. 

6. Physical Processes 

Now we will take a closer look at the physical processes which could cause 
a distortion of the spectrum. Here we will only discuss Compton scattering 
although free-free emission and double-Compton deserve an equally thor- 
ough treatment. 
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6.1. COMPTON SCATTERING 

6.1.1. Collisional Boltzmann Equation 

One may describe the state of the primeval gas of photons and electrons 
in terms of the the density of particles in phase space, i.e. momentum 
and position space. Here we arc not interested in the polarization state 
of electrons and photons so we average over the two polarization states ^ 
It is convenient to measure the phase space density in units of /i = 2-^% 
which gives the the quantum mechanical occupation number, and for 
photons and electrons, respectively. The evolution of can be described 
by the collisional Boltzmann equation which has the form 

^^=0(pj (15) 

where C(p-y) is the scattering term which describes the interactions with 
other particles. Here is a convective derivative along the photon's tra- 
jectory in phase space, while the right-hand-side gives the collision integral. 
If there were no collisions then the Boltzmann equation states that the oc- 
cupation number remains constant along photon trajectories. ^ Included 
in this convective derivative are the all the effects of gravity on the pho- 
tons, which include many of the effects which produce anisotropy. We will 
not discuss these effects further as they are covered extensively in Bunn's 
lectures. 

The collision integral for Compton scattering of unpolarized particles 
after averaging over the polarization state if scattered particles is of the 
form^ 

X (1 - nc(pc)) nc(Pe) (1 + n^iVi)) 
-(1 - ne(Pe)) ne(Pe) (1 + n^{v'-y)) (16) 

Compton scattering in an inhomogeneous medium will produce some polarization 

of the photons, which can be measured, and also effects the anisotropy at the several 
percent level. See Melchiorri and Vittorio this volume. 

^This is true for a phase space defined by a position, x*, and it's canonically conju- 
gate momentum, pi, n(p) measures the particle density with volume measure: d'^'x^dPpi. 
In general relativity there is both the covariant momentum, pi, and contravariant mo- 
mentum, p*. If one measures the density of particles per unit d^x'd^p' the Boltzmann 
equation as expressed above does not apply! 

^This form is determined by the principle of detailed balance which results from the 
time- reversal symmetry of the S-matrix (or classical or quantum mechanics) [45] . 
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where we have (or wih) use the notation 



p^, = ^h p'^ = ^n' |n| = |n'| = 1 E = ^ {m^c^Y + {apef 
Pe = (mec)7/3 (3 = m i = -j2==. (17) 

In eq 16 the values of Pg and e' is determined by energy-momentum conser- 
vation. The 2nd term in square brackets describes the scattering p^, + Pc ^ 
P7 + Pe while the 1st term results from the inverse process, p^y + Pe — *■ 
p^ + Pe. The 1 + n-y factor represents the increased scattering rate due to 
the stimulated emission of the bosonic photons, while the 1 — is the Pauli 
blocking factor giving the exclusion principle for fermionic electrons. The 
factor of 2 in the prefactor counts the two polarization states of the incom- 
ing electrons. The factor c(l — n-(3) in eq 16 is a measure of the relative 
velocity between the ingoing electron and photon.^ Of course, is the 
differential Compton cross-section^ 

Note that this form of the collision integral guarantees that a thermal 
distribution is a fixed point. Substituting a Fermi-Dirac distribution for the 
electrons and a Bose-Einstein distribution for the photon, i.e. 

ME) = --ET^ ^ n^(e) = r-r^ ^ — 7 (18) 

' exp(^ + /xe) + l exp(^+/X7) + l 

will cause the integrand of the collision integral to zero so long as the 
temperature is same for both. Here He, and fj,j are (dimensionless) chemical 
potentials given by the total electron and photon density, each of which is 
conserved by Compton scattering. We expect such a thermal distribution 
to be a stable fixed point since it is the highest entropy state and entropy 
increases according to Boltzmann's H-theorem[45]. In the contexts we are 
interested in the density of electrons is sufficiently low that /Xc ^ 1 and 
Fermi-blocking is unimportant so we may set 1 — ng 1. In this limit 
the equilibrium distribution for the electrons becomes a simple Boltzmann 
distribution, i.e., 

ne(^)=exp(-^-Me) . (19) 

''This relative velocity factor is really determined by the definition of the cross-section. 
The factor is equal to -^Ivi — V2P — ^ '^2^ which reduces to the above expression 

when one of the particles is massless. If both incoming particles are non-relativistic then 
it reduces to the "usual" definition of relative- velocity: |vi — V2|. 

^The outgoing particle momentum pe and p-^ are described by six numbers however 
four are fixed by energy and momentum conservation. The differential cross-section is a 
function of the remaining two parameters, which in this case we have taken to be the 
outgoing photon direction, n'. Any two parameters would do! 
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Henceforth we will ignore Fermi-blocking. 

We are not really interested in the scattered electrons, so we may "in- 
tegrate out" the electron distribution function. The idea is that we know 
the electron distribution function a priori - which is often is a true since 
Coulomb scattering is usually very effective in thermalizing the electron 
momenta. Thus we may rewrite the collision integral as 



cC(p^) = I d% 



s{pLp^) (1 + %(p^)) 



-5(p^,p;)(i + n^(p;)) 



(20) 



where 



'S'(P7,P^) = 



(271^) 



3 I C?''perae(Pe)(l-An)^2^(Pe,P7n) 



c/2 



(21) 

and A(pe,p-y,n') gives the fractional change in the energy determined by 
energy-momentum conservation, i.e. is the solution to the equation 



meC^y + |cpep + cIp^I 



= ^J (meC2)2 + |cpe + Cp^ - 6 (1 + A)n'|2 + c |pe| (1 + A). (22) 



A unique solution always exists with A G [— l,oo). 

We know that the CMBR is very nearly isotropic today, and it is rea- 
sonable to assume that the background radiation was alway isotropic. Since 
we are interested in changes in the spectrum and not anisotropy we may 
also average the collision integral over n to find the mean change in the 
spectrum. Performing the two averages n and n' the collision integral be- 
comes 



C^(e,A) = j dA 



(l + A) 



-S{e,A){l + n^{e{l + A)))nyie) 



where 



5(e,A) = ^^ii±^ [d^n /d2n'5(-n.-(l + A)n'). 

Att J J c c 



■ (23) 



(24) 



To obtain eq 23 we have used a little trick of changing the variable of 
integration for inverse scattering from A to ttx ~ 1' renaming this 
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new dummy variable A. If one looks closely at eq 23 one can see that the 
total photon number is preserved by scattering no matter what the form of 
S{e,A). 

6.1.2. Fokker-Planck Equation 

One important property of cosmological Compton scattering is that, at the 
low redshifts we are interested in, the background radiation photons have 
much lower (total) energy and are moving much faster than the electron 
they are scattering off of. One is bouncing a very light object (the photon) 
off of a much more slowly moving heavy object (the electron) and energy 
and momentum conservation dictates that that energy of the light object 
is nearly unchanged by the scattering (consider bouncing a ping-pong ball 
off of a bowling ball). The electrons are not infinitely massive nor are they 
completely stationary so that the photon energy will change slightly in each 
collision. All this will be reflected in the fact that ^(e, A) when considered 
as a function of A will be a very narrow function sharply peaked around 
A = with width much less than unity. In contrast the A-dependence of 
n-y(e(l + A)) and S{j^,) is a much smoother function in the sense that 
they do not vary much over the region in A where 5'(,A) is significantly 
non-zero. Thus it should be a good approximation to Taylor expand the 
integrand of eq 23 in A about A = 0, but excluding the rapid dependence 
through the 2nd argument of S and truncating at a given order. This is 
a kind of Fokker-Planck equation^. If we expand to 2nd order in A the 
Boltzmann equation becomes a partial differential equation (see eq 8 of 
ref [1]) 



Dn^ 1 d 
where 



e3QeA^^ + (^-A + 2A^+^eA2') (1 + n. 



(25) 



A^ 



iVe(7T J-1 



dAA"5(e,A) A"' 



NeaT J-1 



de 



(26) 



Fokker and Planck actually considered the case where the momentum is only slightly 
changed in each scattering and proposed Taylor expanding to 2nd order in the small 
change in momentum. For Compton scattering the direction of the photon will change 
significantly so the momentum change is not small, but the energy change is, and ex- 
panding in the small fractional energy change. A, is an obvious generalization. It is useful 
to consider expanding to higher order than 2nd. 
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and we have used the electron density, A^e; introduced the Thomson cross- 
section^, (TT, and defined the Thomson optical depth, r: 

(27) 

This optical depth gives the expected number of Compton scatterings of 
low energy photons off of non-relativistic electrons. 

The form of the equation is reminiscent of a diffusion equation which 
is good description of the physics, the small changes in photon energy at 
each scattering causes the photons to diffuse in energy space. The 
term causes a net drift toward increasing energies while the (1 + n^)nj 
will cause a net drift to lower energies (if it's coefficient is positive). We 
expect these drifts and diffusion to sum to zero in thermal equilibrium, 
i.e. when n-y has a Bosc-Einstcin distribution (eq 18), the electrons have a 
Boltzmann distribution (eq 19), and the two share a common temperature. 
This consideration alone suggest that for a thermal electron distribution 
with temperature that we should expect 

-2l + 4A^ + eZ^ e ^^^^ 



A2 

Another feature of eq 25 is the differential operator in front, which 
guarantees conservation of photon number. This will persist to all order in 
the A-expansion. In fact one can pretty much guess the 2nd order Fokker- 
Planck equation without knowing much about the Compton cross-section. 
We will take a more constructive approach below. 

6.1.3. Compton Cross-Section 

To compute the Compton collision integral, or it's Fokker-Planck approxi- 
mations one needs to use the Compton cross-section. The relativistic (Klein- 
Nishina) differential Compton cross-section in an arbitrary rest-frame is[l] 

dV 3(7T 1 - P'^ 



d^n' 167r [i _ H'-p + aj-^{l - n-n')]^ 

(1 - J-)(l - fi-n') 



1+ 1 

{l-h.p){l-h'-(3)J 



(1 - h'-p) [1 - h'-f3 + a7-i(l - n-n')], 

'^Thomson scattering is the non-relativistic and classical limit of Compton scattering 
as first described by J.J. Thomson. 
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Figure 8. Plotted is the distribution of fractional energy changes, A, experienced by low 
energy photons scattering off of an isotropic distribution of electrons with velocity /3c. 
The left panel shows the distribution for (3 = 0.01, 0.02, and 0.05; while the right panel 

shows the distribution for 13 = 0.1, 0.2, and 0.5. For graphical clarity wo have adjusted 
the heights of the curves to have unit amplitude at A = 0. The maximum and minimum 
values for A are dictated by energy and momentum conservation. The distribution is 
narrow and symmetric for small /3 and becomes wider and more skew for larger velocity 
electrons. This positive skewness gives the heating of the photons by the electrons. The 
Fokker-Planck equation approximates the photon distribution function by the first few 
terms in it's Taylor series about A = when convolving with these distributions. This 
is liable to be a good approximation for scattering off of low velocity electrons since the 
A-distribution is sharply peaked around A = 0. 



where 



For many astrophysical applications, and especially those related the the 
CMBR there are two small numbers which enter this cross-section. Firstly a 
is very small for the microwave photons wc see observe today, roughly 10~^. 
Clearly as we go to higher redshifts the background photons become more 
energetic, but a remains small in the redshift range relevant to the CMBR 
spectrum z < 10^. The 2nd small number is /3 since we are almost always 
interested in non-relativistic electrons. If one is interested in a thermal 
electron velocity distribution then a small /3 expansion is equivalent to 
a small expansion. In most applications the a ^ so we will 

concentrate more on the higher order terms in and not a. 

To proceed it is probably easiest to follow the methodology of Barbosa 
[1], where one expands the cross-section in a but not p. For a thermal elec- 
tron distribution one can compute the moments, A", in the Fokker-Planck 
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expansion analytically, and only at the end one should Taylor expand the 
result in Te about Te = 0. One may rewrite eq 24 as 



(31) 



where ^ gives the differential cross-section wrt to the fractional change in 
photon energy. So for example, expanding everything to zeroth order in a 
(which we denote by the superscript ^°^) we find 



(l-/3-n) 



dA A 



= aTi^(A,/3sgn(A)) 



(32) 



where 



F(A,6)=sgn(A)x7^(l-i^^) 

"3(1 - b- f{3 - l>^){2 + A) ^ (1 ^ b){l + A) 



1666 

3(l-fe^)(2fe-(l-6)A) 
3266(1 + A) 



In 



1 + 6 
(4(3 - 36^ + 6^) 
+2(6 + 6 - 66^ - 6^ + 26^)A 
+(1-62)(1 + 6)A2) , (33) 



and H(J is the Lorentz-Heaviside function which is unity for positive argu- 
ment and zero otherwise. We see that this function is only non-zero for 



AG 



2/3 2/3 



1 + /3' 1-/3 



(34) 



and, as promised, for small /3 is sharply peaked around A = 0. We plot this 
function for various values of /3 in fig 8. 



6.1.4. Moments of A 

With this general expression for ^^^^ given above one can compute, to 0th 
order in a, the A-moments which are the coefficients in the Fokker-Planck 
equation (some of these may be found in ref [1] : 



^=(0) 4 



= 4 



10 



meC" 



+ 



( k% 
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=(0) 2 



A2 = ^4^2(5 + 16/32) 



mcC^ 



A3™=^7^/?4(21 + 23/32) 



+ 47 
252 f kT^\'^ 



O 



+ 



=(0) A M / kT \ ^ 

= -J-78/34(i47 + 1554^2 + gSQ/J^) = ^ ^ ] + O 

525 5 V'TT'eC 



meC^ 
fcTe 



(35) 



and we also find that A"^' 

=(0) 



(0) 



since 



has no dependence on e. The 



fact that A'^ = 1 tells us that, to 0th order in a and all orders in j3 the 
scattering rate per unit volume is cNgaT-^ The coefficients in the Fokker- 
Planck equations are determined by the average of the electron velocities 
indicated, and these expressions hold whether or not the electrons are in 
thermal equilibrium. For a thermal distribution these velocity moments can 
be computed exactly in terms of modified Bessel functions [1], however we 
have found it convenient to expand these functions to the appropriate order 
in temperature. It seems that a Taylor series to a given order in A is less 
accurate than the same order Taylor series expansion in Tg. To keep track 
of the various terms in the expansion let us devise the notation 



0(n, m) = 



(36) 



There are no terms ~ 0(0, 0). One finds that 



=(m) 



:(m) 



A2"-i ~ A2» ~ 0{n, m) 



(37) 



so to include all the terms of order ~ 0{n, m) in one must make a Fokker- 
Planck expansion to order 2n in A. It is probably not worthwhile to go 
to high order in these expansions, since one can circumvent this expansion 
by doing the collision integral. Nevertheless the first few terms give useful 
analytical expressions. 



*The total (Klein-Nishina) cross-section starts to fall below the Thomson cross-section 
when the center-of-mass photon energy rises to close to meC^, i.e. when 7a > 1. A careful 
look at eq 29 will show that by setting a = in this equation we are ignoring terms of 
order 07. For microwave photons this approximation should be good for computing the 
total cross-section as long as 7 < 10^ i.e. for anything less energetic than ~ 500TeV elec- 
trons. In contrast to compute the small effects on the spectrum from Compton scattering 
one should include 1st order terms in a whenever a > (3^. 
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6.1.5. The Kompaneets Equation and Relativistic Corrections 

The lowest order non-zero Fokker-Planck equation, given by the expansion 

of eq 25, is the Kompaneets equation 



dn.y 



]_d_ 



de 



+ 



{1 + n. 



0(1,0) 



0(0,1) 



(38) 



where the order of the two terms are indicated. This equation was first 
pubhshed by Kompaneets [42] in 1957 and probably developed earlier as 
part of the Soviet thermonuclear weapons program. For hotter gas one can 
add terms 0{2, 0) which yields an extended Kompaneets equation[72] 




nieC 



(39) 



Further terms in this expansion will be derived in ref [72] although it is not 
clear how useful they will be since extensive numerical work has been done 
with the more accurate collision integral (e.g. ref [56]). 

6.1.6. The Generalized Sunyaev-Zel'dovich Effect 

The idea of the Sunyaev Zel'dovich distortion is that one starts out with a 
background radiation which is close to a blackbody spectrum, just what we 
expect to be produced by the early universe, and it is slightly distorted by 
the action of hot ionized gas through the Compton scattering process we 
have just described. In this small distortion limit we need just substitute 
in a blackbody spectrum, n^^ of eq 2, into the right-hand-side of the Kom- 
paneets equation. Let us generalize this idea a bit by instead considering 
the more general Fokker-Planck expansion which is an expansion in and 
a. In this small distortion limit the different terms will add linearly to the 
total distortion which we may write as a sum 



\n^ 



EE 

n>Om>0 



X = 



(40) 



where 



"V/"(n,m 



dT 



kT, 

rUeC^ 



kT. 



7 



mgC^ 



(41) 
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Figure 9. Plotted is the small deviation in intensity from a blackbody divided by the 
classical y-parameter caused when blackbody photons pass through a hot gas of electrons. 
This is computed using the extension of the y-distortion given in the text. The gray band 
is centered on the classical j/-distortion which applies when kTc <C m^c? . The black lines 
are for electron temperatures of 1, 2, 5, 10, 15, 20, and 25 keV. We have of course assumed 
To 3> T^. The curves intersect at the zeros of the function Ansz°\ 



and the superscript correspond to the 0{n,m) contributions to the 
Fokker-Planck expansion. Substituting n^^{x) into the various terms of 
eq 39 we find that 

A</'(x) = 

^ ' (e^ - 1)2 \^ 2 - 1 5 (e^ - 1)2 

7 3 (e" + l)(e2- + 10e" + l) ^^ 
+ 10^ (e--l)3 ]• 
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One does expect that to each order in energy that a blackbody spectrum 
is a stable solution when the electron and photon temperature are equal so 
one should expect the sum rule 

N 

J2^nir-"\x) = (43) 

n=0 

and this does seem to be true for iV = and and N = 1. 

The classical Sunyaev-Zel'dovich y-distortion contains only the 0(1,0) 
and O{0, 1) terms and may be written 

XC^ / g-^ -|- ]^ \ 



where 



J rrieC^ 

This is the y-distortion plotted in fig 5 and used in eqs 11&12. To sec how 
much this classical formula errs we compare the different expression for a 
range of electron temperature in fig 9. We see that the 0(2,0) corrections 
become significant when Tg > 5keV. This 2nd order distortion agrees well 
with the computation of the collision integral by Rephaeli[56] so higher 
order corrections do not seem to be important for Tg < 15keV. 



7. The Future 

In the past decade we have witnessed astounding advances in the measure- 
ment of the CMBR spectrum. After decades of tantalizing evidence of devi- 
ations from a blackbody spectrum we find that the spectrum is amazingly 
close to a perfect blackbody. No longer is it possible to consider a universe 
with a very hot inter-galactic medium, or that hydrodynamic forces could 
have played a large role in forming the large scale structure. There is also 
little room for non-equilibrium energetic events in the early universe at red- 
shifts < lO'^. In a way this is most unfortunate. The thermal equilibrium 
state contains the least information - all remnant of events in the universe 
before z ~ 10'' have been thermalized to nothing, or more precisely to one 
number: the temperature. At the moment we really don't know how to in- 
terpret this number, other than to make a rough comparison to the number 
of baryons which is observationally rather ill-determined. Perhaps some day 
we will have cosmogenic theories which will predict the baryon-to-photon 
ratio with great accuracy. 

Observationally we are approaching a brick wall which is the Galaxy. At 
the present level of sensitivity Galactic contamination from dust and syn- 
chrotron radiation is an important contaminant at all wavelengths. Galaxy 
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modeling which makes use of a spectral and spatial structure of the Galaxy 
observed at a variety of wavelengths will improve as sensitivities improve 
however there will be a limit to how accurately one can subtract off the 
Galaxy even given perfect data. We won't be making observations outside 
of the Galactic plane any time soon! 

Yet there is still a lot of room for improvement on the decimeter and 
meter scale anisotropies. Also there is this tantalizing evidence for negative 
spectral deviations 

Things are not bleak. In fact spectral distortions of the CMBR is a 
rapidly growing field. Multi-frequency observations is beginning to be the 
norm for CMBR anisotropy experiments, and with the CMBR satellites we 
can expect literally millions of measurements of the CMBR spectrum in 
different directions on the sky. Admittedly there is a big difference between 
absolute measurements of the CMBR flux and differential measurements 
which are required for anisotropy since the anisotropy spectral measure- 
ments are modulo any DC spectral distortion. However it is just his sort 
of measurement which will make improved Galaxy subtraction possible. 
The spectral information obtained will tell us mostly about the Galaxy 
and extra-Galactic radio sources, however with millions of measurements 
one can always hope for something a little more interesting. Along these 
lines there is the cluster S-Z effect which is a rapidly growing field. With 
increased sensitivity we can look forward to S-Z selected cluster catalogs, 
measurements of radial cluster velocities through the kinematic S-Z effect, 
and these studies can work their way down to galaxy groups and even large 
scale structure filaments of hot gas. We can even hope to measure the gas 
temperature from spectrum if it is hot enough. In the future we can ex- 
pect the spectrum and anisotropy measurements to become increasingly 
intertwined. 
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